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Abstract 

We investigated the light propagation by means of the Robertson-McVittie 
solution which is considered to be the spacetime around the gravitating body 
embedded in the FLRW (Friedmann-Lemaitre- Robertson- Walker) background 
metric. We concentrated on the time delay and derived the correction terms 
with respect to the Shapiro's formula. To relate with the actual observation 
and its reduction process, we also took account of the time transformations; 
coordinate time to proper one, and conversely, proper time to coordinate 
one. We applied these results to the problem of increase of astronomical 
unit reported by Krasinsky and Brumberg (2004). However, we found the 
influence of the cosmological expansion on the light propagation does not 
give an explanation of observed value, ciAU/cit = 15 ± 4 [m/century] in the 
framework of Robertson-McVittie metric. 
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1. Introduction 

The appearance of radar /laser gauging techniques and spacecraft ranging 
has enabled us to measure the distance from the Earth to other inner planets 
and the moon directly and precisely. Recent years, these techniques have 
been drastically improved, for instance, the planetary radar measurement 
achieves the observational accuracy of the interplanetary distance within a 
few 100 [m] , the spacecraft ranging a few [m] and the lunar laser ranging a few 
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[cm]. In these observations, the round-trip time of hght/signal is measured 
by the atomic clocks on the Earth. Nowadays, the advent of laser cooled 
atomic clocks has led to the improvement of realization of SI second and it 
is expected that this advance will go on. Therefore the arrival time measure- 
ment is the most accurate observation among the all kinds of astronomical 
and astrophysical ones. 

Due to the improvement of these techniques, it is important to develop 
the rigorous light propagation model which contains the various general rela- 
tivistic and relating effects, and this subject is actively investigated by many 
authors, e.g. 0, M, H fl H H H M, 0, H Q and references therein. 



These theoretical developments also play a crucial role in testing the gravi- 



tational theories 43, 44 



On the other hand, an arrival time measurements have also devoted to the 
improvement of lunar and planetary ephemerides, such as DE [42] , EPM 33 



VSOP [2I, and INPOP 11], and the determination of astronomical constants. 
Especially among these constants, the astronomical unit (of length) AU is 
one of the fundamental and important one which gives the relation of length 
units; [AU] of astronomical system of units and [m] of SI ones. AU is currently 
determined from the arrival time measurement data of planetary radar and 
spacecraft ranging within the accu racy of 0.1 [m] or 12 digits level as 1 [AU] = 



1.495978706960 x 10^^ ± 0.1 [m] [33|. Therefore from the point of view of 



the fundamental astronomy, the development of a accurate light propagation 
model is of great importance. 

Until now, the light propagation formulae have been derived in terms of 
the post-Newtonian approximation, g^j^iy = rjfj,^ + h^^, \h^u\ -C 1 in which the 
background metric r]^jj is the static Minkowski one. However, it is interesting 
to examine the influence of cosmological expansion on the light /signal prop- 
agation since the accuracy of astronomical observation is rapidly growing. 
Of course, it is presently hard to observe such a influence in the solar system 
experiments. Nevertheless, the progress of these or relating measurement 
techniques may enable to detect its trace in the future date. 

In this paper, we will study the light propagation passing near the massive 
star as the Sun which is embedded in the cosmological expanding background. 
We will adopt the Robertson-McVittie spacetime introduced in section [2], and 
concentrate on the gravitational time delay and derive the extra delay due to 
the effect of cosmological expansion with respect to the standard formula by 
Shapiro in section [31 Further in order to relate with the actual observation 
and its reduction process, we will take account of the time transformations 
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in section IH from coordinate time to proper one in subsection 14.11 and con- 
versely, from proper time to coordinate one in subsection 14. 2[ As the ap- 
plication of results, we will consider the problem of increase of astronomical 
unit recently reported by 25| in section [51 



2. Robertson-McVittie Spacetime 



The standard form of Robertson-McVittie metric is given by |35l . |27| . Il4l . 



GM 

2c^ra(t) \ ^2^^2 



1 _|_ GM 



+ 1 



GM ~ 
2c^ra{t) 



a'{t){dr' + r^de^ + r^sm^ 



where G is the gravitational constant, M is the mass of central body, c is 
the speed of light in vacuum, and a(t) is a scale factor. ([1]) is regarded 
as the spacetime around the point mass singularity embedded in the FLRW 
background metric, and coincides thoroughly with the Schwarzschild solution 
in the isotropic coordinate when a{t) = 1, and the FLRW model for curvature 
parameter k = when M = 0. 

However, recent dynamical model of planetary motion (EIH equation of 
motion), the various light propagation models cited in the previous section 
and other observational models in the solar system are formulated not in the 
comoving coordinate system as ([T]) but in the barycentric celestial reference 
system (BCRS) or corresponding to BCRS based on the post-Newtonian 
framework jioj. And the various astronomical constants are also derived 
in BCRS or corresponding to BCRS; for example the secular increase of 
AU, dAU/dt that we discuss in section [5] is evaluated based on barycentric 
dynamical time (TDB) Therefore in order to discuss these effects and 
cosmological one in the same framework as far as possible, it seems to be 
more adequate to take the reference system as close as possible to BCRS. 
Hence we adopt the following radial transformation and convert ([1]) into the 
nearly proper coordinate system 35 , 14 , 15 , [l^, 

3 HQ, 



R = a{t)r 1 



GM ~ 
2c^ra{t) 



(2) 
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We notice that the matching of BCRS and cosmological reference system is 
important problem in the fundamental astronomy and this issue is discussed 
by (l9|,[2lj. Using ([2]), (P is rewritten as, 



, ^ 2GM\ 2 , 2 



1 _ 2GM 



here H{t) = a{t)/a{t) is a Hubble parameter and the overdot denotes the time 
derivative. If H{t) does not change, namely, the Hubble constant at present, 
H{t) ^ Ho = (/io/3.08) X lO^i^ [1/s], ho = 0.7, then we can introduce the 



time transformation to remove dtdR term in ([3]) 35, 14 



„ Hn f RdR , , . 

cT = ct + — / 7 — , , (4) 



2GM 



In a practical sense, the Hubble parameter H{t) varies with time. Nonethe- 
less, when focusing on the solar system experiments, the actual round-trip 
time of light/signal and the time interval in which the observational data 
is stored are much shorter, at most t ~ 100 [yr], than the age of Universe, 
Tu ~ 10"^ [yr]. Thus, let us assume H(t) changes adiabatically as. 



(t - to). (5) 

It is suited to put, dH/dt\Q ~ Ho/Tu ~ 10^^^ " yr)], hence, it follows, 

X 100 [yr] ~ 10-^^ [1/s]. (6) 



Ho ~ W-'' > ^ 
dt 



From ([H]), Hq produces the dominant effect of cosmological expansion. Sub- 
stituting (jl]) into ([3]) and limiting the equatorial motion 9 = Tr/2, we obtain, 

,2 2GM Hi 2\ 2,2 A 2GM H^ ^\ ^ ^ 2,^2 

ds^ = -\ l ^r^ c^dt^ + 1 + + -^r^ dr^ + r^dcfP, 7) 

y c r c J \ c r c ) 

where we expanded the coefficient of dr^ and replaced T —>■ t and R r, 
respectively. 
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Here we note that the assumption H{t) = const, leads Robertson-McVittie 
solution to Schwarzschild-de Sitter model, see e.g. 19|, |35j. Friedmann equa- 
tion is, 



2 kc^ 



8nG A 



(8) 



where k is the curvature parameter, p is the density of universe, and A is the 
cosmological constant. When the cosmological background is the de-Sitter 
universe, p = and k 



0, we have. 



constant, 



(9) 



in this case a{t) = exp(cyA/3t). Replacing Hq by iQ and inserting ([2]) and 
dl]), into ([T]), we can recover the Schwarzschild-de Sitter solution. 

As a consequence, although our metric ([7]) is equivalent to Schwarzschild- 
de Sitter model, we will begin with this metric as the first attempt of our 
investigation. 



3. Time Delay in Robertson-McVittie Spacetime 

The world line of light/signal is the null geodesic, ds"^ = 0, then from ([7]) 
it results, 



ct{r, ro) 



dr' 



1 + 



2GM 



+ 



,.!2 



\1 



X 



2GM _ Hi ^,2 

(P-t' (P- 

_ 2GM 
c^r' 



Hi ^,2 



2GM 



Hi ^2 



-1/2 



(10) 



in which ro is the closest point (impact parameter) between the light/signal 
and the central body. Remaining 0{GM / ,Hq/ (?) terms in the integrand, 
we obtain. 



t(r, ro) 
trpN{r, ro) 



osmo 

1 



(11) 



rl 
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tcosmo(r,ro) = gy'r2-rg(2r2-r2), (13) 

here tiPN is the Shapiro time delay in the 1st post-Newtonian approximation, 
and tcosmo is the extra one caused by the cosmological expansion. If we 
assume the Earth (E) and the objective planet /spacecraft (R) are almost at 
rest during the round-trip of light / signal, then the net of round-trip time T 
becomes, 

T = TiPN + Tcosmo = 2[t(rs, ro) + t(r^, ro)]. (14) 
The time delay produced by the cosmological expansion is. 



Cosmo 



3c3 



2 r 



- ri (2r| - rl) + ^rj, - rg (2r| - r^) 



(15) 



4. Time Transformations 



4.I. Transformation from Coordinate Time to Proper Time 

(1141) and fll5p give the time delay in the coordinate time. However, the 

measurement of round-trip time is actually carried out by the atomic clocks 

on the Earth which ticks the proper time r. Therefore, we must convert f|T^ 

and f|T5|) into the observer frame. 

To this end, it is sufficient to consider the equation of proper time in the 

quasi-Newtonian approximation, 

^ = l-^-^-Mr^ (16) 
dt 2c2 cV 2c2 ■ ^ ' 

Taking the orbital radius and the velocity of the Earth, rE and ve-, respec- 
tively, the measured round-trip time T becomes. 



T — Tlpn + Tcosmo- (17) 



E 



Making use of (fT^ , the part due to the cosmological expansion relating with 
i^o is, 



rcosmo = ^{^[v'4^(2r|-roV\/4^(24-ro^ 



rl-rl + 



rR-ri 



.— ^r||21n 



te + 



r p - r, 



^0 




+ 2 M^o 



+ 



r2 (2r 



(18) 



This is the extra time delay observed by the atomic clocks on the Earth. 
0{c~^) terms in first and second line are the leading ones due to the cosmo- 
logical expansion and 0{c^^) ones from third to sixth line are the coupling 
terms i.e. cosmological expansion and Shapiro delay, and so on. Hence the 
coupling terms are considered to be next order of dominant terms in terms 
of post-Newtonian approximation. 



4-2. Transformation from Proper Time to Coordinate Time 

Next, we consider the transformation from proper time to coordinate 
time. This transformation is also of importance when obtaining the position 
and velocity of planet / spacecraft (they are the reflectors of light / signal) re- 
ferring the proper time r of observer on the Earth (practically, coordinated 
universal time, UTC or international atomic time, TAI). This is due to the 
fact in the ephemeris the position and velocity of celestial body is calculated 
as the function of the coordinate time t which is the independent variable of 
the equation of motion. 

Applying the Keplerian energy integral. 



GM 



CM 
2ao ' 



(19) 



where is the semi-major axis of planetary orbit, and (1161) is rewritten as, 
assuming the Earth (observer) moves along the circular orbit, r = rE = ao; 



dr _ 3GM 
'dt~ ~ 2cVij 



(20) 
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By integrating fl20l) we have, 



T'lPN 



osmo; 

(21) 

-^)(^-^o), (22) 



TCosmo — 2(^^^~^^^' ^^"^^ 

here Tipn represents the relation between the coordinate time t and the proper 
time r, sometimes called the time ephemeris (the accurate analytical treat- 



ment of TiPN is given by, e.g. [28|, l29|, l9j). rcosmo is caused due to the 
cosmological expansion. 

5. Application to Increase of Astronomical Unit 

In this section, let us apply above results to the increase of astronomical 
unit recently reported by Krasinsky and Brumberg (2004) jisj. They found 
that from the analysis of high precision planetary radar and martian space- 
craft ranging data (the arrival time measurement), the astronomical unit of 
length AU increases with respect to meters as dAU/dt = 15±4 [m/century]. 
This reported values is about 100 times larger than the present determina- 



tion error of AU, see 33| and section [T] in this paper. The similar variation 



of AU is also corroborated by [42 



Before entering the AU issue, we briefly summarize how the AU is derived 
from the analysis of observational data. These days, the AU is determined 
based mainly on the planetary radar and the spacecraft ranging data since 
these data give the interplanetary distance £obs directly and precisely as, 
^obs = ctobs [m] where tobs [s] is the observed round-trip time of light/signal. 
While, the planetary ephemerides provide the interplanetary distance rtheo 
in the unit of [AU] as the theoretical value. Then to compare rtheo in [AU] 
with observed tobs in [s], Ttheo is converted into ttheo in [s] as, 

AU 

^theo = '^theo— = rtheoTA [s] (24) 

in which = AU/c is called the light time. Practically this step looks around 
the various effects, e.g. the relativity, the solar corona, and the troposphere, 
see 5.32 of |37|. Therefore the light propagation model plays a key role 
in deriving the AU. We also note that rtheo depends on a constellation of 
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astronomical constants and parameters via the equation of motion and some 
other relations. Then AU is optimized and derived with other parameters 
simultaneously by the least square method, as satisfying tobs = ^the 
However, when Krasinsky and Brumberg rewrote ttheo as. 



tieo ■ 



^ (au . fir) . (25) 

and fit to observational data, they discovered dAU/dt is non-zero and positive 
value dAU/dt = 15 [m/century], where T is the time interval counted off 
from some initial epoch. Notice that estimated dAU/dt does not mean the 
expansion of planetary orbit and/or the increasing of orbital period of planet. 
According to Krasinsky ji^], the observations do not show any traces of 
such kinds. Further the determination error of inner planetary orbits in the 
latest planetary ephemerides is also smaller than the observed dAU/dt = 15 



[m/century], see e.g. Table 4 of [33j. Therefore the observed dAU/dt may 
relate with not the dynamical aspect of planetary motion but the light/signal 
propagation. 

Therefore, let us examine if the extra time delay derived in previous 
sections can explain the observed dAU/dt. We first apply (fT8|) and evaluate 
the Earth-Mars ranging concretely. We suppose the Earth (E) and the Mars 
(R) move along the circular orbit of radius, te ~ 1.5 x 10^^ [m] (1 [AU]) 
and r/j ~ 2.28 x 10^^ [m] (1.52 [AU]), respectively. Figure [T] shows the 
estimated extra time delay Tcosmo in the proper time of observer. The figure 
is plotted as the function of relative angle between the Mars (pn and the 
Earth (pE, \4>R ~ 0b| for \(f)R — (j)E\ > vr/2. The order of 'Tcosmo is at most 
about 3.2 X 10^^^ [s] which is considered as the typical magnitude of 'Tcosmo 
for the planetary radar and the spacecraft ranging. Here let us evaluated the 
order of magnitude of each part in f[T5]) : the leading terms H^r^/c^ ~ 10~^*, 
the coupling terms GMHqt'^/c^ ~ v'^Hlr^/e' ~ 10~^^, and H^r^ /e' ~ 10"^^ 
where we took tejVe of Earth as r,v. However, it is much smaller than the 
current measurement limit, that is, the internal error of atomic clocks, 10~^ 
[s]. 

Next, we calculate (1231) presuming again the Earth moving along the 
circular orbit, te = 1.5 x 10^^ [m]. The difference between the coordinate 
time and the proper one caused by the cosmological expansion is, in the 
interval of 100 [yr], 

rcosmo = -^^(^ - ^o) ~ -4.1 X 10-'° [s]. (26) 
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Figure 1: Extra time delay Tcosmo measured in the proper time of observer on the Earth. 



This is also much smaller than the relative error of atomic clocks, 10 ^ [s]. 
Tcosmo causes the apparent orbital variations of longitude of planet, but it 
is S(j) = noTcosmo ~ -8.2 X 10-27 [rad] where uq = ^GM/r% ~ 2.0 x 10~^ 
[rad/s] is the mean motion of Earth. Then rcosmo does not has much effect 
on the observed dAU/dt. Rounding up the results, the effect of cosmological 
expansion on the light propagation dose not give the explanation of observed 
increase of AU in the current solar system experiments. 

Here we mention the attempts of Krasinsky and Brumberg (2004). They 
also investigated the secular increase of AU in terms of cosmological expan- 
sion deriving the approximate metric of Einstein equation, (5) of [25| which is 
the modified version of the standard form of Robertson-McVittie metric ([T]). 
In such a comoving form, they examined both the planetary motion and 
light propagation, then found there exist the large orbital variations in the 
radial and longitudinal directions due to the cosmological expansion ((19) 
and (20) of [25|], respectively). However, these dynamical perturbations are 
completely canceled out by the time transformation between the coordinate 
time and proper one relating with light/signal propagation ((22) and (23) of 



^(5) of [25|] can be recovered from the lowest order of H]) with respect to GM/c?, 
replacing the time coordinate dt a{t)dt and a{t) y/~A, see appendix A of [25l |. 
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25|, respectively). Hence they concluded that the cosmological effect cannot 



account for the observed dAU/dt. 

Before closing this section, we shortly note the dynamical effect due to the 
cosmological terms in the planetary motion. This problem has been studied 
by many authors 0, H 0, 0, fl, 0, S, H Ei, fl, S 0, [ii, 0, Q • And 
they showed that the cosmological influence on the planetary dynamics, e.g. 
the additional perihelion advance, is negligible small. We also confirmed that 
the cosmological terms are too small to give an account of the increasing of 
AU. In appendix [Aj we briefly summarize the results of the cosmological 
effect on the planetary motion based on ([7]). 



6. Summary 

In this paper, we investigated the light propagation passing near the 
gravitating body as the Sun embedded in the cosmological spacetime by 
means of the Robertson-McVittie solution. We concentrated on the time 
delay and derived the correction terms with respect to the Shapiro's formula. 
We also took account of the time transformations; coordinate time to proper 
one, and conversely, proper time to coordinate one. These treatments are of 
importance when dealing with the observed round-trip time and its reduction 
process. 

As the application of results, we considered the problem of increase of 
astronomical unit reported by Krasinsky and Brumberg (2004), since the AU 
is currently derived from the analysis of the round-trip time of light / signal. 
However, the effect of the cosmological expansion is too small to produce the 
observed dAU/dt = 15 [m/century] and then it dose not give the explanation 
of increase of AU. 

Here it is worth to mention the interpretation of dAU/dt. Although we 
and Krasinsky and Brumberg (2004) investigated the secular increasing of 
astronomical unit in terms of the cosmological effect, another possibility is 
suggested, that is, increasing of AU is arisen due to the lack of the calibrations 
of internal delays of the radio signals within the spacecrafts and this may be 
the most plausible reason of observed dAU/dt. Nonetheless, till now, the 
origin of dAU / dt is far from clear therefore this issue should be explored by 
means of the every possibility we can imagine. And the re- analysis of dAU/dt 
using new data sets is also expected. 

As we mentioned in the end of section[5l so far, the dynamical terms due to 
cosmological expansion in the planetary dynamics has been studied actively. 
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while the influence on hght/signal propagation has hardly been examined. 
However, the accuracy of astronomical/astrophysical observations, especially 
the optical, radio and arrival time ones, rapidly increases in the solar system 
and they may make it necessary to incorporate the cosmological effect in the 
future date. For this purpose, our investigation in present paper is the flrst 
attempt to this direction. To consider this problem, the matching of BCRS 



and cosmological reference system discussed by |19| . |2l[ is also important 
issue to be investigated. 

([7]) expresses the completely static gravitational fleld, and is equivalent to 
Schwarzschild-de Sitter model. Nevertheless, as long as we focus on the solar 
system experiments, our results may be regarded as the dominant effects 
due to the cosmological expansion, see in However in order to study 
the cosmological influence on the light / signal propagation in more detail and 
general cases, we must treat H{t) and/or a{t) as a function of time and adopt 
the time dependent spacetime model. This subject will be investigated at 
some future time. 

We would like to express our gratitude to the referee for fruitful com- 
ments and suggestions. We acknowledge to G. A. Krasinsky for providing 
the information and comments about the AU issue. We also appreciate T. 
Fukushima, M. Kasai, H. Asada, Y. Itoh and M. Takada for fruitful discus- 
sions and comments. This work was partially supported by the Ministry of 
Education, Science, Sports and Culture, Grant-in-Aid, 18740165. 

A. Planetary Motion Based on Metric ([Tj) 

Cosmological influence on the planetary motion has studied by many 
authors as we mentioned in section [51 Then, in this appendix we limit to 
summarize our results based on ([7]) briefly. 

The equation of motion in the quasi- Newtonian form becomes. 

Supposing initially circular orbit r = rg = const., the mean motion is ex- 
pressed as, 

n = ^:^no + 6n, 5n = -r^o^ ~ 1.0 x IQ-^^ [^^^/^j^ (28) 
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where no = yGM/rQ ~ 2.0 x 10 [rad/s] and we took the orbital radius of 
Earth te as tq. From fl28|) the variation of longitude (p is evaluated in 100 
[yr] as, 

S(f) = 6n{t - to) 3.5 X 10"^^ [rad]. (29) 

Also using fl28l) . the orbital period of planet (in this case that of Earth) is 
written as, 

T = ^c^To + ST, 5T = ~ 7 X 10-2^ [s], (30) 

here Tq = 27i /rio. Finally putting r = + 5r and inserting into fl27j) . it 
yields, 

^4o-2 / ^3 7:72 \ 

- W - ^) ^ H. (31) 

Therefore we find that the estimated orbital variations 6r, 6(j) and change 
of orbital period 6T are much smaller than the observed dAU/dt = 15 
[m/century]. 
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